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The problem of electronic structure of random metallic 
alloys with its varied applications and extensions has attracted 

the attention of condensed matter theoreticians. The foundation 

of successful theory was laid by Soven <:i967> and Kirkpatrick et 
al C1970> through the Coherent Potential Approximation CCPA> 
Soven crystallised the idea of single site CPA, with the 

intention of introducing complex site energies Cor potentials> 
which characterised the effective medium in .a self- consistent 

manner. Within multiple scattering formalism this meant that if 

a single exact potential is embedded in this medium, there is no 
extra scattering on the average. 

The simplest case f or which most of the CPA calculations 

have been done is that of binary random alloys without 


short- ranged 

order and clustering : 

A B 

c i-e 

where c 

is 

the 

concentration 

of A atoms In the 

alloy. 

Usually, 

the 

site 

energies are 

taken to be random Cdiagonal 

dlsorder> 

and 

the 

overlap Integrals are assumed to be 

non-random, 

the 

band 


widths are reasonably equal. Also it Is assumed that correlated 
scattering from two or more sites is ignored. This is known as 
single site CPA or ICPA. 



V 


This scheme has been ext-ensively applied in bobh bhe 
bi^hb- binding and bhe Koni^inga-Kohn-Rosboken <KKR> scabbenin^ 
approaches. Ib has been increasingly felb bhab bhe single sibe 
CPA appears bo be Inadequabe in describing bhe ePPecb ojf bhe 
clusbersj> exbended disorder and shorb-ranged order. 

A new formalism is proposed Cin chapber II> Por developing a 
selb-consisbenb CPA encompassing bobh single pobenbials as well 
as clusbers of pobenbials wibh perhaps shorb ranged order. This 
mebhodology is an exbension bo bhe random sysbems of bhe 
embedding bechnique inbroduced by Inglesfield <1977^ 1978, 1981> 

for impuribies and surfaces. This involves bhe parbibion of bhe 
Hilberb space on which bhe Hamilbonian is defined inbo bwo 
subspaces, I and II, wibh an inberface S bebween bhem. The 
advanbage here is bhab bhe bime independenb Schrodinger equabion 
in bhe f uU Hilberb space may be solved in bhe smaller subspace 
I. The effecb of bhe resb of bhe subspace, t.e., Cbhe complemenb 
of subspace I > is incorporabed in bhe effecbive Hamilbonian on 
subspace I as a surface pobenbial. This surface pobenbial can be 


consbr ucbed f rom 

bhe Green 

funcbion 

of 

bhe subspace II 

which 

sabisfies 

bhe usual boundary 

condibions 

ab 

infinity 

along 

wibh 

Neumann 

bovindary 

condlblons 

on bhe 

inberf ace S. 

This 

Green 

f tmcbion 

depends 

upon bhe pobenbial 

in 

the 

subspace 

II. 

Here, 


subspace I is described by random binary alloy A B , whereas 

e 1 — c 

subspace II is described by an effecbive non-random pobenbial, bo 
be debar mined self-consisbenbly applying CPA or KKR-GPA depending 
upon bhe sysbem bo be sbudied. 



In chapt-er III t^his: embedding CPA f ormalism las applied *to 

"the Coherent Jellium models;. Here subspace II is replaced by 
translatlonally symmetric^ energy dependent^ site independent 
Coherent Jellium medium. The following problems are considered : 

I. Electrons with s-symmetry^ in square well potentiaal. 

II. Electrons with d-symmetry> in square well potential. 

III. Electrons with s-symmetry in -V/r potential, 

IV. Electrons in non- spherical potential well. 

In chapter IV one dimensional Embedding- KKR formalism is 

developed to calculate the electron density of states of one 

dimensional array of flat potentials within a single well CPA and 
a double well CPA. Here the Gx*een function of subspace II can be 
calculated f rom the one dimensional analog of KXR-CPA theory 

developed by Butler <1976 > to study the electronic properties of 
transition metals and their alloys. 

In concluding chapter further extensions and usefulness of 
the embedding formalism have been discussed. 
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CHAPTER I 


INTRODUCTION 

Ever since *the enunciat^ion of" -the Bloch’s -theorem C1928> 
considerable successful efforb has ^one inbo bhe understanding of 
electronic properties of crystalline^ periodic systems. 
Concerted and extensive effort began into the study of 
compositionally disordered crystalline solids^ i.e.* i-s^ndom 
alloys > only in the late sixties. The kind of problem one 
encounters in random systems are also multifarious. For example 

I. The understanding of various disorder linked 

order-disorder transitions, magnetic transitions etc. 

II. The theoretical understanding of various electronic 

properties like density of states, optical properties, 
photo- emission, conductivity, dielectric response and 
electron related magnetic properties. 

III. The understanding and description of the nature of 

disorder-driven localization transitions of electrons 
in random systems. Nature of localised states and 
transmission in disordered networks. 

IV. The understanding of the elementary excitations in 

random systems such as phonons, magnons, polar ons 


etc. 



In t.his -thesis we sh«all ^ddx^ess our^selves exclusively to the 
undei^standing: of the density of elect3>onic states of random 

metallic alloy systems. 

As we have experienced in the case of ordered systems^ once 
the Hamiltonian model is set up> the solution of the quantum 
mechanics is essentially the solution of the Schrodinger 
equation. However> with the introduction of disorder, we face a 
new problem, i.e., the lack of sufficient information about the 
sample in the Hamiltonian. The potentials which describe a 
disordered solid are characterised by random parameters: random 
in space as in quenched disordered solids or random in time as in 
thermally disordered systems, f or example, electron in contact 
with phonon bath. Dirty alloys Cthe so called Mooij systems > at 
hi^h temperature may involve both spatial and temporal disorders. 
A particular realisation of these parauneters, either in a g^lven 
sample or at an instant is what we shall call a ^confil^uration^of 
the system. Properties of a given ^configuration^ not be 

pertinent. We often take recourse to configuration averaging. 

Before we go into the study of various methods of 
configuration averaging, it is important to examine why 
configuration average is needed and what is to be averaged. 

To understand the first part, i.e.» why it is needed, let us 
consider the specific case of a crystalline substltutlonally 
random binary alloy, in which the regul ar lattice sites are 
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x^andomly occupied by A oz^ B bypes of" aboms:: randoTn GuNi Is a ^ood 
example. At. zero bemperabure bhe randomness is quenched and may 
be described by a seb of" random occupabion variables n. , which is 
debermined by baking bhe values O and 1, whebher a sibe labelled 
by i is occupied by an A or a B bype of* abom. A configurabion is 
debermined by a parbicular assignmenb o£ O's and l^s bo bhe 
sequence n^. A given sample corresponds bo a given configurabion 
whereas differenb samples bo dif"f*erenb conTigurabions. There are 
2^ differenb conrigurablons for samples conbaining N aboms In a 
binary alloy. This is a very large number for macroscopic 
sysbems. Vhen an experimenbaJisb balks aboub bhe characberisbics 
of a sample he measures, he is obviously nob inberesbed in bhe 
variabions bebween bhe large number of configurabions. Rabher he 
seeks average brands among his samples. If bhe flucbuabions of 
a characberisbic from configurabion bo configurabion is 
negligible compared bo bhe mean, bhen ib is bhe configurabion 
average which should be compared wibh experimenbally debermined 
characberisbics. 

The above discussion illusbrabes why configurabion averaging 
is imporbanb. Ib also indicabes bhab ib is bhe physically 
measured proper by which should be configurabion averaged: bhe 
densiby of sbabes or bhe response f uncbions rabher bhan bhe 
Hamilbonian or self-energy, bhe diffusion probabilibies rabher 
bhan wave ftincbions. Anderson C1969> gives a lucid argumenb on 
bhis bopic bo explain some of bhe anomalous resulbs on bhe nabure 
of bhe specbrum of random bighb-blnding Hamilbonians, which were 
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being pi?opos;ed ab bhab bime phased on configtii^abion averaging. 
Becausse of bhis^ bhe Green function technique is far more 

suitable for bhe study of disordered systems than the 

Schrodinger eqpoation approach. 

The above description also clarifies that configuration 
averaging is a valid procedure only when the probability 
distribution of the property under consideration is sufficiently 
well behaved so that the average dominates over the higher 
moments, otherwise it is not profitable to talk in terms of the 
configuration averaging. A classic example is the intensity of 
star light transmitted through layered media with randomly 
varying ref ractive Indices studied by Chandrasekhar <1960!> and 
the closely related problem of the resistance of a disordered 
chain. Kumar and Coworker <1985,1986 > have shown that the 
variance of the resistance of a disordered chain diverges much 

faster than its mean as the length of the chain increases. It is 
not meaningful therefore to talk about the average resistaunce of 
a very long disordered chain. There are indications that the 
same statement may be valid f or three dimensional disordered 

systems, certainly on the localised regime and probably even when 
the Fermi energy lies Just above the mobility edge. 

Now it^s worth-while to say something about spatial 
ergodicity. For a large enough system, in the limit of Infinite 
size, all possible environments are achieved even in a single 
sample, as shown in figure <1.1>. 




Figure 1.1 : Representation of different 

environment within a sample. For n 

•all environments are realised. 
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In a global- pi^obe sample all possible envix'onmenbs in bhe 
sample are weigrhbed wibh appropriate possibilities. A global 
property is thijs automatically conTiguration averaged^ even if it 
refers to a single sample. This sampling of the entire 

configuration by a spatially global probe is analogous to the 
sampling in statistical mechanics of all configuration space by a 
dynamically evolving system. Hence the name spatial ergodicity. 
Such global properties need not be configuration averaged all 
over again. This idea is implicitly behind a large body of the 
"‘large cluster" calculations involving either isolated or 
immersed clusters. One example is the application of the 

Recursion method to large disordered clusters CHaydock, 1987>. 

There are various mean field approaches to study the 
conf iguration averaged one and two particle Green functions f or 
disordered systems^ which are directly related to the averaged 
density of states and response functions. But during past two 
decades Coherent Potential Approximation CGPA> has been one of 
the main theoretical development in the area of study of 
electronic structure of random alloys. The simplest case for 
which most of the CPA calculations have been done is, that of a 
random binary alloy. Here the site energies are taken to be 
random Cdiagonal disorder> and the overlap integrals are assumed 
to be non-random, i.e., there is no off-diagonal disorder. Also, 
it is further assumed that the correlated scattering from two or 
more sites is Ignored. This is known as single site CPA, 
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The v^dldlt^y of t^he single sife CPA -theory has; been 

recognised through application to various model systems as well 
as to realistic model Hamiltonians suitable for disordered 
metallic alloys. The need to j^o beyond IGPA has been felt 

because of the following drawbacks of single site CPA : 

I. It takes exactly into account only single*-slte 

scattering and correlated scattering from clusters 
is not accounted for> 

II. The effect of potential fluctuation is automatically 

suppressed while doing such an average, 

III. It cannot properly take into account the 

off-diagonal disorder in the Hamiltonian, and 

IV. The vertex correction Ci.e., back-scattering effects>‘ 

vanishes in the single band case as an artifact of 
short-ranged potentials in single site CPA. 

The most obvious extension of the theory to overcome these 
limitations is to go beyond the single-site approximation by 
doing a general formalism so that the effects of scattering from 
statistically coupled clusters of various sites and that of the 
off-diagonal disorder could be included. Nickel and Krumhansl 
C1971> tried to make such generalisation using a method based on 
the correlated cumultant expansion to find an electronic density 



of sfafes of a one dimensional model of an alloy. 
BuflerClP72,1973> also t-X'ied fo calculate bot-h fhe i>eal and 
ima^inaT'y panfs of bhe Oreen function using: a self-consisfenb 
cluster mefhod and bobh of bhese mebhods a 3 >e inbended bo imp3>ove 
bhe CPA by ^oing beyond bhe single- si be approximabion bhnougrh 
making applicabion bo one- dimensional big^hb binding models. Such 
a ^eneralisabion nob only posed severe analytic problems such as 
having: disconbinuibies ab various energies bub also led bo 
mulbi-valued density of states. 

It was pointed out by Haydock etal <1972>> Mookerjee C1973>^ 
Muller- Hartman <1973>, Mills St Ratanavaraksa <1978> that for any 
real potential function^ bhe Green function must have certain 
mabhemabical properties, called Herg^lobz properties, A complex 
function fCz) is defined as Herglobz if 

i> Im fC 2 > < 0 for Im z < 0 
il> Singularities of fCz> lie on bhe real axis, 
lii> f< 2 > — > l/E as E — co Cwhere z ■ E + 10^> 

In order bo g^eb physical results, Herglotziciby shotild be 
related in any approximation bo bhe Green f unction. The 

fundamental difficulties with most theories that were proposed 
during bhe late 1960^s and early 1970*s as extensions of bhe CPA, 
are that they lead bo approximate Green f unctions which are nob 
Herglobz. lb is now realised that bhe preservation of bhe 
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t^he development. oJf a t^heory which includes scai.t.er*ing from 
clusters > shorf -ranged order > off-diagonal disorder^ or fhe 
effects of positional disorder. 

There have been attempts to include some of the ef f ects of 
scattering from clusters of atoms^ at least conceptually by using 
a technique called the Molecular Coherent Potential Approximation 
CMGPA> as discussed by Tsukada C1969> and Butler C1973>. This 
method is a straight-forward generalisation of the CPA approach 
in which the lattice sites of an alloy are partitioned into cells 
that contain more than one lattice sites. Though in this method 
the Herglotz analytic property is preserved, the main objection 
for MCPA not being a satisfactory theory is the imposition of 
fictitious unit cell boundaries, and consequently, translational 
invariance is lost. 

All the above techniques described do not prove a rellad:>le 
solution to the theoretical problem of developing a 
self- consistent theory of the electronic states of an alloy that 
goes beyond the single-site approximation. The problem of going 
beyond the single site CPA theory has been proposed by Mills 
<1978> in his original studies Cknown as Travelling Cluster 
Approximation > on the problem of a Herglotz average Green 
f unction within the f rame work of diagrammatic perturbation 
theory. But it^s very cumbersome to evaluate the contribution of 
the higher order terms f or a real three dimensional lattice by 


Mookerjee C1973 1975 a^b>c > inbx-oduced a ^enei^al 

formalism known as Au^menbed Space Formalism CASF> bo calculabe 
configrurabion avera^ginig of a general funcblon of random variables 
by usin^ absbracb operabor f ormalism of qtianbum mechanics. He 
was able bo esbablish bhe idea of incorporabing scabbering 
effecbs from sbabisbically coupled sibes in a finibe size 

clusber, as well as bhe Imporbance of off-diagonal disorder in a 

« 

self-consisbenb manner. This formalism came oub wibh success in 
producing Herglobz analyblc Green f uncbion bhrough appUcabions 
bo various model sysbems. 

On bhe basis of bhe Parbibion bheorem and bhe Embedding 
ideas of Inglesfield we develop a mebhodology bo sbudy bhe 
elecbronic sbrucbure of random sysbems. The equabions for bhe 
Green funcblon are confined bo a finibe clusber region and bhe 
influence of bhe resb of bhe region appears as a surf ace 
pobenbial. The Herglobz analybic proper bies are sabisfied in 
bhis case. We illusbrabe bhe mebhodology by applying ib bo a 
binary disbribubion of dlf f erenb single wells in a Coherenb 
Jellium. This is bhe alloy generalisabion of bhe impuriby work 
of Inglesfield referred bo in bhe bexb. This formalism is 
developed in mosb general f orm bo bake Inbo accounb clusber 
effecbs as well. Exbended defecbs and impuriby clusbering in a 
random alloy can also be baken care of wibh ease. This 
mebhodology is a very good bechnicpae bo deal wibh liquid alloys 
wibh shorb-ranged ordering. 



Her»e we also develop a one dlmens:ional analog or t-hls 
Embedding jTonnalism for single as well as pair po-ten-tial binary 
alloy impuribies in a one dimensional array of flab pobenbial 


wells. 



CHAPTER II 


AN EMBEDDING FORMAUSM TO THE COHERENT POTENTIAL APPROXIMATION 


Ono or -the most. powei>rul techniques availedble dunin^ the 
last two decades f or the study of random metallic alloys is 
Coherent Potential Approximation CGPA>. This scheme has been 
extensively applied both in the tight-binding and 

Korringa-Kohn-Rostoker CKKR> scattering approaches. But it has 
been extensively f elt that the single site Cor single potential > 


CPA appears to 

be inadequate 

in describing 

the effects 

of 

clusters. 

extended 

disorder 

and 

short- ranged 

order which 

may 

dominate 

situations like 

order- disorder 

transitions. 

for 

instance. 

In this 

chapter we 

wish 

to propose a 

new formalism 

f or 


developing a self-consistent CPA encompassing single potential 
wells as well as clusters of wells and short-ranged order. This 
methodology is an extension to random systems of the Embedding 
technique introduced by Inglesfield C 1971,1972, 1981 > for 

impurities and surfaces. Before going to develop this formalism 
let us very briefly review Coherent Potential Approximation. 


2.1 : THE COHERENT POTENTIAL APPROXIMATION 


The electronic density of states per atom per spin for any 
system is directly related to the single particle Green function 
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nCE> ■ - i_ im Tr GCE+iO*> 

TIN 

C2.1.1> 

where n is *the t^ot-al nuwnber of a'toms in t^he sys-tem. The Green 
opera’tor Cor resolvent, operat.or> G is ^iven by bhe Hamilbonian H 


GCZ> * CZI-H> 


C2.1.2> 


where Z is a complex variable and I is t,he idenbiby mabrix. 


For any disordered sysbem bhe elecbronic density of states 
will be ^iven by trace of the configuration averaged single 
particle Green operator <G>. If the random variation of the 
potential is similar throughout the system, i.e., randomness is 
homogeneous, then the configiuiration averaged electronic density 
of states will be 

nCE> - - — Im Ti><GCr,E+10'^» 

TTK 

C2.1.3> 


where r is any site vector in the system and the averaged Green 
function will be 

<G> m 


C2.1.4> 

The CPA is an attempt at an approximate evaluation of this 
configuration averaged Green function. It maintains the analytic 
features of the exact Green function and interpolates correctly 
between several individual limiting cases, the Virtual 

Crystal Approximation, atomic and dilute limits, weak and strong 


disorder limits. 
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There atre four different^ approaches t-o €he CPA as described 
in Ihe expensive literature on the approximation . 

I. Diagrammatic approaches using propagator and locator 

formalisms introduced by Yonezawa and Matsubara C1966>, 
Yonezawa C1968>, and Leath C1968,1970>. 

II. The effective medium approach introduced by Seven 

a967>. 

III. Multiple scattering ideas^ first introduced by Anderson 

and McMillan C1973>, Mookerjee <197S> and Kaplan &t al 

C1980>. 

IV. The Augmented space approach by Mookerjee <1973> and 

Kaplan et al C1980>. 

All different approaches yield the same final result within 
the single site approximation. This indicates that the single 
site CPA is unique. In this thesis only the effective medium 
approximation will be described. 

2.1.1 : THE EFFECTIVE MEDIUM APPROXIMATION 

The idea of the effective medium approximation <Soven 1967> 
gives a simple physical insight into the CPA. Here we replace 
the actual random potentials at different sites by an effective. 



energy dependent, translation^dly symmetric potential; such that 
the Green iTunctlon corresponding to the effective Hamiltonian is 
equal to the configuration averaged Green function of the actual 
system under consideration. By this definition the effective 
Hamiltonian is related to the configuration averaged Green 

operator as 

<GCZ» - <ZI-H 

C 2.1.5 > 

From equation C2.1.4> and C2.1.5> it is clear that the 
constructed effective medium is described by complex energy 
dependent potentials sitting at various atomic sites. Though the 
effective medium does not correspond to a real physical system, 
its Green function is related to the configuration averaged 
properties of the random system. 

Let us illustrate the procedure in a simplified single- band 
tight-binding model with diagonal disorder. Our Hamiltonian in 
the tight ““binding basis is 


H « 


E 

.tv 


V 




t J 


tj 


T. . 

tj 


C 2.1.6 > 


where 1, j are the site indices, and are projection [ iXi | 

and transf er | i>< j | operators respectively on the space spanned 
by the basis < j l». The diagonal or site energies are random. 
Ve take V. . to be some averaged nonrandom value <V. >, 


Now we define oux' effective Hamilt^onian in €he seune 


‘tigrh't-bindlni^ basis as 


H _ « r Z CE> P. r <V > T. 

.V t VJ IJ 

^ j 


C2.1.7a> 


In a single sibe CPA bhe self energy Z^CE> is in general 
complex and energy dependents. Ib is also non random and 
independenb of bhe sibe index» i.e., Z^CE> » Z^CE>. 


Wibh bhis our efjfecbive Hamilbonian becomes 


H , «z CE> r p + r <V. > T. 

*f’f O V . . l-J 

V V J 


C2.1.7b> 


The only unknown cjuanbiby in bhe calculabion of bhe elecbron 
charge densiby is Z^CE>. The problem is now bo debermine bhis 
self energy or effecbive medium energy self-consisbenbly^ which 
yields 


G CZ> * < G<Z> > ■ CZI^H 
eff eff 




To do bhis we consider an exacb pobenbial^ which is embedded 

at. the site within the effective medium tPleur*e 2.1> In such 

a way that no extra scattering is produced on the average. 
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Figure 2.1 tSingle potential embedded m 
effective medium. 
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In mat.hemat.ic^Ll language 


<0'^> - G , - < G > 


< 2.1.9a > 


g'^ « C zi-h" ) ‘ 


where > 

«X1 

< 2.1.9b > 

is -the Green function for bhe medium in which all bhe slbes are 

occupied by Z^<E> except the site I which is occupied by s , . The 

corresponding Hamiltonian <h'' > is 

*1 f 


H „ = H „ + <£ -Z > P 

C2.1.10> 

Substituting equation C2.1.10!> in equation C2.1.9b>> it can be 
written as 


G*^ = r ZI-H , - Cs. - E >P 1 
[_ I O i J 


-1 


and svibstituting equation C2.1.8>, we obtain 


-i 


- °.rf [*- ] 

<:2.i.ii> 

Expanding t-he right, hand side of t.he eqiiat^ion C2.1.11) we get. 

0^-0 ,^+G fc^.-Z >p1g , Ji-G >p.11 


-I O ^ G 


C^.-Z >P. fl-G , -z >F 1 G 

1. O V »ff t o vj «ff 


C2 i 
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Subst-it-u-ting 


b - <£ -Z >P, f 1-G {c£ -Z >P 11 

I V O t V o Ij J 


-1 


in equation C2.1.12> we obtain 


G + G t G 

•If ©ff i off 


^nd t^aking conri^\ira€ion average 


<G > m o + O <t > G 

eff erfl I &ii 


In order t^o sa-tisfy equat^ion <2.1.9a> we writhe 


<t > = O 

V 


C2.1.13:> 


C2.1.14> 


Equat^ion C2.1.14> ensures t.hat. when an extra potential is 
embedded in the effective medium there is no extra scattering^ on 
the average where t^ describes the scattering:- 

The self energy Z^CE> can be obtained from this equation, 
solving it self- consistently. Here we are Interested in studying 
the electronic properties of binary alloys^ for which the 
site either can be occupied by an A atom or a B atom. Therefore, 
for binary alloys without short-ranged order, equation C2.1.14> 


reduces to 



c 


+ <1“G> 


O 


O 


1-<G><^ - Z > 

A O 


B 


1“< GXiC > 
B O 


<:2.i.i5> 


where c is t.he concentration of A atoms in the alloy. This' is 
the CPA equation for random binary alloy. To solve this equation 
seLT-consistently a suitable initial choice or is made for 

<G> to st30?t with. Ve then calculate Z^CE!> f rom the equation 
C2.1.15> and then from Z <E> we calculate a new G. CZ>. This goes 

O V 

iteratively until Z^CE> converges. Care must be taken to 

preserve the Herglotz analytic properties CMooker jee,1973> 
related to the physical necessity that the energy spectrum be 
real and the density of states be positive. This must be 
preserved at all levels of iterations. 


2.2 : A METHOD OF EMBEDDING 


The method of Embedding was introduced by Inglesfield C1981> 
to study the electronic properties of impurities and surfaces. 
This methodology involves the partition of the space on which the 
Hamiltonian is defined, into two subspaces I and 11 with an 
interface S between them. This Is shown In the figureC2.2>. 


The Hamiltonian for the full space is 



2 ; 


where Vj ^nd U are •the potentials in the subspaces I and II 
respectively. 

If ^<r> is the wave function defined for the full space then 
the Schrodinger equation of the system will be 


H ^Cr> * E ^Cr> 

C2.2.1> 

where E is the energy eigenvalue. 

The advantage here is that the above time independent 

Schrodinger equation in the f uH space may be reduced to a 

Schrddinger equation in the smaller subspace I, through an 
effective Hamiltonian in subspace I. The effect of the rest 

of the subspace II, i.e,, the complement of I, is incorporated in 

this effective Hamiltonian as a surface potential, V ,as shown in 

B 

figure C2.3>. 

The effective Hamiltonian Cdefined in subspace I> is. 


- T+V,+V 
1 s 

< 2 . 2 . 2 > 


This surf ace 

potential can 

be 

constructed from 

the 

Green 

function G Cr,r^> 

of siibspace 

II 

which satisfies 

the 

usual 


boundary condition at infinity alon^ with Neumaiui boundary 
condition on the surface the normal derivative of the 


Green function vanishes on the surface S : 


2 


“hn 


O Cr^rO I 
o » 




fi 


S 


O 


C2.2.3> 


To solve ’the Schrodingei> equation in notion I we find the 
ejipectation value of the energy of the wave function ^Cr> of the 
full space> defined in subspace I as a trial wave function ^Cr> 
and in subspace II as the exact solution of the Schrodinger 
equation at some energy with the condition that the solutions 
match correctly on the boundary surface S: 


+U<r> >V'Cr> * 0> r in subspace II 


C2.2.4:> 


and 


yyCr > « <pCr > 

& St 

C2.2.5> 

The expectation value of the Hamiltonian of the full 
subspace is 


Jd®r Cr) 


< 2 . 2 . 6 > 


Siibst^it-ut-ine equation C2.2.4> in equation C2.2.6) we obtsdn 


f dr^y/* Cr>v'Ci'> 

I II 


E ■ 


i fdV p*<r >f ^ 1 

* ■’ * • I- a; — s:: j 


dn 


r d^r p Cr>joCr>+r d r Cr>v^Cr> 



The firxal surface integral in t.he numerator comes from the 
discontinuity in the derivative between <^Cr> and ^<r> on the 
surface S. We now express the normal derivative across 

the boundary surface S in terms of <^Cr by using the Green 

St 

function G Cr^O for region II, satisfying the following 

o 

equat.ion 


C- i +UCi'>-jc> G Cr-^yO ■ <5Cr-i«0 

2 r o 


where r and r^ are in region II. 


C2.2.8> 


Multiplying equation <2.2.8> by v^CrO and equation C2.2.4> by 
G Cr,r’> and subtracting one from the other we get 


C-i V^u;Cr»G Cr.r^^ +C- V^G Cr,rO>V'<»'> ■ -<5{:r-r'^>V'<r> 

2 r o 2 r- o 

C2.2.S» 

Int,ecratlne equat^lon C2.2.9) through the region II yields 


vCr> - - r d®r"EG Cr,r07^,V/Cr0-\y€r0'7*,G Cr^OI 

^ 2 o r o 


This integral can be evaluated using Green theorem 


J 


If Cr>V^ gCr>-“gCr>V^ f Cr>-ld^r 




fCrD 


dgCr> 

dn 


gCr> 


df Cr>‘J 
dn 

a -* 



C2^.10> 



2 ‘ 


wlier*e Is t*he normal deriva'Live aL 'bhe surface S direcLed 

oui>wards from inside t^Iie volume. Usln^ Green t^lieorem t^lie volume 
in’be^ral of equaf-ion C2.2.10D reduces f-o trlie following^ surface 
integral 


-i f |g Cr^rO - y/<rOdG ] 

^ \ ‘ » -^T- “ J 

B dn 

8 

C2.2.11> 

If we construct G^ to have zero derivative on S, this gives 
us an equation relating the amplitude of the wave function in II 
and its normal derivative on S : 


V/Cr> - --f G Cr ,r«0 du/CrO 

2*' 9 0 a 

9 


dn' 


snd putting r on S 


> ■ --f d^n' G Cr ,r'> dv^Cr*) 

s 2J s o s a a 


dn' 


C2.2.12> 

The inverse of this gives dv^Cr >/dn in terms of V'Cr which 


a a 


we have set equal to ^Cr >, as given by the equation C2.2.13> 

9 


dv/Cr > ■ - af d^r' ^Cr'> KCr ,r'> 

9 9 S 8 9 


dn 


C2.2.13> 


where the functional inverse K of G is defined by 

O 


f KCr,r"OG Cr^»,rOd^r" 

J o 


5<r-rO 
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Now t-o reduce the volume integral J d®r 

to a surface integral, we vary the energy £ in the equation 
C2.2.4!!>, but still ensuring that v/Cr> matches with ^Cr>, over the 
in'terf'fiice S, "to ob'tain 

S£ v'<r> + £ 6v^Ci'> 

Multiplying^ this ;and the complex conju^gaite of" equation <2.2.4> 
r-espectively by y/ and 6v^Cr*> and subtracting^ one f^rom the 

other and integrating’ through subspace II ^ the above volume 
integral yields; 

C2.2.15> 

The volume integral on the right hand side can be reduced to 
the surface integral using Green theorem as follows 



1 

2 



Cr > 


<r > 

fit 

“Si 


* 

Cl' > 
8 


dn 


dv'Cr > 
8 

&£ 



C2.2.16> 


Substituting dy/Cr >/dn and it's complex conjugate and using 

6 S 

the symmetry of the Green f unction in the surface integral 
C2.2.16>, the volume integral becomes 


r d r 4* Cr>4‘Cr> ■ -fdr rdr'<^Cr>-x KCr ,r'> ^Cr'> 
J J a J B a a£ m a a 

11 a a 


C2.2.17> 
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Stjl>s;'blt#ut.in^ ©cjuat-ion C2.2.13> ajnd i&qiaait^ion C2.2.17^ In equ^'Lioii 
« C2.2.7> we obt^adn 


E » 


f d^r ^*Cr>H<(>Cr>+ - fd^r p^Cr > 

^ 2 J B a 

s 


dn 


^ J J <^CrO 


^ fdr rdr*^^Cr>~ KCr <5i><r^> 

J aJ B ^ B &£ B B ^ B 

B B 


r* 3 ^ 

I d r ^ Cr‘>^Cn> 

•^i 

-f d^r r dV’ <p* ex' > ^ KCr ,r»> <?!>Cr'> 
a B as & B B 

S B 


C2.2.18> 

This expr^ession gives E puGrely in berms oT a “trial wave 
Tunebion <^Cr> debined in bhe subspace I and in bhe inberJTace S. 
Now bhe variational minimizabion of" E with respect bo bhe trial 
wave Tunebion (p Cr>, i.e.^ SEySp Cr>»0, gives bhe effective 
Schrodinger equation. If we evaluate bhe Green f unction at bhe 
eigenvalues s ^ E on bhe spectrum, this effective Schrodinger 
eqfuabion reduces bo 


-- + - 6Cn-n > 

2 r- 2 a 


^ + v<:r>^r5 

B 


+ SCn-n >r dV' KCr ,E> ^frCrO - E 

s J s a s a 

s 

wl-th r In subspace 1. 

C2.2.19> 

This effective Schrodineer equation acts only on subspace I 
and the effect of the subspace II appears as a surface potential 



operator 6Cn-n >r dV' KCr The term Cl/^><5<:n-n ><.&yan > 

B ** 9 m 9 9 

9 

ensures t.ha*t €he effective Hamilt-onian remains Hermit^ian in 
region I alone. Usually rei^ion I is a small subspace of" bhe 
Hilbert, space^ so bhat equabion <2.2.19> is a signiricanb 
simplificabion> provided bhab we can debermine bhe surf" ace 
pobenbial on bhe surface wibh ease. 

2.2.1: THE CALCULATIONAL PROCEDURE 


For acbual calculabions of bhe resolvanb operabor^we need a 
represenbabion of bhe effecbive Hamilbonian Cin subspace I> in a 
suibably chosen basis. To obbain bhis we expand bhe brial wave 
funcbion in subspace I in berms of counbable basis funcbions 

X Cr> which span subspace I and reflecb bhe symmebry of bhe 

m 

solubion we wish bo sbudy 

<^<rl> * Z a ;;t 

m mm 


C2.2.20> 


Using: bhis basis and minimizing energy E wibh respecb bo bhe 

var iabional pao^ameber a , we obbain bhe mabrix represenbabion of 

m 

the effective Schrddinger equation C2.2.19>: 

ZH a«EZO a 

m rvm m m nm m 

C2.2.21> 


The effective Hamiltonian matrix is 


H 

nm 


jy 


r X 

n 




V* +VCr>1 


X + 

2 


; 1-*“ 


V Cr > ^x ^ 
m ^m m 


an 



28 




X y K<r 

ns s 


r' E> X Cr* 

, a, 8 


and the overlap matrix Is 


<2-2.22a> 


O 

nw 




X 

n 


X <r'> 

m 


C2.2.22b> 


2.3: CPA THROUGH EMBEDDING 

We shall now develop CPA 'thorough this Embedding technique. 
The effective medium is described by an effective unknown 
potential term. This can be constructed from the Green function 
G^<r where U is the coherent potential of subspace II ^ to 
be determined self- consistently. 

Bef ore going on to develop the CPA through embedding let us 
review the philosophy which underlines the CPA as well as the 
Embedding method. 

Inglesfield^s idea is based on the fact that the solutions 
of the second order elliptic differential equations Csuch as 
Schrodlnger equation> in all space may be solved by confining 
ourselves to a smadl finite closed subspace. The effect of the 
remaining space is incorporated in the effective Hamiltonian as a 


surface potential. 
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In the single site CPA we find the self energy or effective 
medium energy self -consistently by embedding an exact potential 
at the site within the effective medium in such a way that 

no extra scattering is produced on the average. In mathematical 
language this is 


<G^<V, 


,U» 


<Vi> 


G 


C2.3.1> 


where G is *the Green Tundion oT t^he errect^lve medium and 

oil 

'the averaged Green tunction in eqxiation <2.3.1> is averaged over 
different configurations of Vj in subspace I alone. 


Preserving the essence of these approaches we now develop 
our new , Embedding- CPA formalism in the following steps: 


We divide the space into two subspaces I and II with an 

interface S. In subspace I we take the potential VCr>«=Vj<r> 

exactly as a random potential whereas in subspace II the 

potential UCr> is replaced by an effective energy dependent 

non-random^ translationally symmetric coherent potential 

U ^Cr,E>: U ^Cr4^R,E> = U , Cr,E>, where R is lattice 

coh coh coh 

translation vector. 


As described in CPA, coherent potential or self energy is 
obtained, embedding a small region of the exact potential with 

a particular configuration in this effective medium Cregton JJ> 
so that no extra scattering is produced on the average : 



C2.3.2> 


<GC)r^%V,U ,E» ■ G Ct>^*,U ,E> 

off 

where »V,U »E5 is 'the embedtiing^ Green f'unction of 

subspace I sabisfyin^ bhe effecbive Schrbdlnger equation: 

CH*-EI> GCr^' ,V,U ,E> - 6Cr-rO, 
is the corresponding effective Hamiltonian. 

The effective medium Green function G Cr^*,U,E> is 
obtained from a Hamiltonian where the Coherent Potential UCr,E> 
is substituted for the entire space. Equation C2.3.2> reflects 
CPA in a more general way without considering a particular basis 
in region I. 

2.3.1: CALCULATION OF EMBEDDED GREEN FUNCTION 


First let us construct the embedded Green function which is 
needed to study the electronic properties of a disordered system 
applying Embedding-CPA method. This Green function satisfies the 
following equation: 


-- 7*QCr^%U ,E> + ^Cn-n > d GCr^%U, E> 
^ " to 


+ C VCr>-E > GCr,r^U ,E> 


+6Cn-n >r dV^» KCr .E> G<r" ,r%U.E> - 6Cr-rO 

a m a B a 

where r^^are in subspace I. 


<2.3.3> 



Let us take a repr^esentation of equation <2.3.3) in a basis 
which spans subspace I and which is suitably chosen to 
reflect any symmetry of the solution required. Expanding the 
Green function in terms of this set of basis functions: 

G<r,r' ) - ]] G x<r’:> ;t:<r' ) 

Ttft\ 

rim 

and substituting in equation <2.3.3) we obtain 

G<r,r>) - E A:^Cr)+ i 6<n-n^) ^ + V<r);t^<r) 

nm L 

a 

+ 6<n-n )r dV" K<r ^ ,E>x <r" )-E;t <r')U <r' ) » <5<r-rO. 

s J s B a coh n s in 

s J 

<2.3.4> 

Multiplying this by ):j^<r););^<r') and Integi'ating over r and r^ 
through the region I we obtain 


Z 

nm 


<H. -E 
kri 


kn 


°nm °ml 


'kl 


<2.3.5) 


where t^he Hamiltonian and overlap matrix elements are ^iven by 
equation C2.2.22a> and equation C2.2.22b> respectively. Linear 
independence of" basis then yields 


E<H -EO 
, nk 
k 


nk 


> O 


km 


6 

nm 


C2.3.6> 


So the matrix representation of" the embedded Green function is 
the following resolvant matrix : 

G 

nm 


<H-EO>" ‘ 
nm 
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In cont^inutim of* en^p^y It# is xisejf ul t#o wopk wit^li 

t#li© local donsit-y of" st#3t#os nCp,E)^ t#hat# is, chtap^o donsit^y 

of elections with enep^y E. This can be dipecbly found from t#he 
ima^inapy pant of t#he embedded Gpeen function: 

nCp,E> * — E I ^ CpD, p in subspace I 

nm L nm J n m 

C2.3.7a> 

The bobal densiby of st#at#es is pelabed bo bhe above bhpou^h 



pc:p>*= ^ J dE nCp,E> 

"CX) 

C2.3.7c> 

where p is bhe chemical pobenbial. The number densiby of 
elecbrons is given by 


n» J dE nCE> . 

-OD 

<:2.3.7d> 

2.3.2: CALCULATION OF THE EFFECTIVE MEDIUM GREEN FUNCTION 


We shall iilusbrabe our formalism by applying it# t#o a simple 
example. Here region I is a sphere of radius r^ cenbred ab bhe 
origin wibh a random single muffin -bln pobenbial V^Cx^y. The 


boundary is a spherical inberface S> as shown in bhe figureC2.4>. 



Figure 2.2 


: Partition of space into two subspaces 



Figure 2-3 ; Advantage of embedding formalism. 



Figure 2.4 


: Muffin-tin potential embedded in region 
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For simplicity of calculation, we take a Coherent Jellium in 
region II. This assumes that the Coherent potential is 

independent of r but strongly energy dependent : UCr,E>« UCE> 


The Green function Q^<r,r',E,U> can be expanded as 


G Cr,r%U,E> 

O 


-i« E 
L 


lC«r JfCxr > . ^ 

I < I a h^Cxr^> 

h^Cxr > 

L a 




C2.3.8> 

where are spherical Bessel T unct^ions, h^ are Hankel 

funct^ions and are spherical harmonics and is defined as 

1/ 2 ^ 

C2<U<E>*~ E > 1 . Here L and r denot^e t>he angular momen-ta unif 

vector in the direction of r, described by the an^ixlar variables 
<l,m> and CO^<p> respectively. 


On the boundary surface S this Green function becomes 


Q Cr ,U^E> 

os s 


-i« E WCj^^h^> Y^Cr>Y^CrO 

hfC«r > 

I S 


where 


E 

L 


g Cxr >Y, <r>Y, Cr' 

I S L. la 


C2.3.5» 


[ h C^’Xjy 
I s 


C2.3.10> 



and t,he Wronskian is 


WCJ ,h ) ■ j Cxr >h 'Cxr > - j /Cxr >h Cxr > 

«.<U s I fil< 


This is rela'ted t,o it,s functional inverse as 


r k Cxr >e <«r > ■ 1 

el e l s 


C2.3.11) 


Substltutinc equation C2.3.10> in equation C2.3.11> we get 


k, Cxr > ■ i i 

I S 4 


X r^ VCJ^,h^> 


h^<>er' > 
I s 

h C«X’ > 

I 5 • 


and 


r-. ys. 


KCr ,U,E>* £ k C«r >Y, Cr>Y, Cr>> 

» s I s L L 


<2.3.12> 


C2.3.13> 


Here t.he « dependence of fhe surf ace pof enf ial shows fhaf 
the effective Hamiltonian or embedded Green function depends upon 
the coherent potential of subspace II. 



CHAPTER III 


COHERENT JELLIUM MODELS 


In i^hlst cKapt.er we shell illus'tnat.e otxr Embedding- CPA 
foiMTialism by simple examples. Ve bake examples of pandom blnapy 
alloys, ^ t.h& concenbpablon of A bype of aboms 

in bhe alloy. binary alloy pobenbial in region I. 

In region II, we bake a Coherenb Jeilium pobenbial, i.e., a 
pobenbial bhab is independenb of bhe spabial variable r, bub 
depends only on energy E UCr,E> * UCE>. In bhis model bhe 

labbice effecbs are suppressed. Such models are appropriabe in 
cases where long range labbice ordering is absenb, as in liquid 
alloys. 


3-1: ELECTRONS WITH S-SYMMETRY IN SQUARE WELL POTENTIAL 

We shall sbudy bhe s wave local denslby of sbabes in flab 
pobenblals, i.e., bhe pobenblals VjCr!> * or in region I 

and zero oubside. In bhis subspace I we choose our basis 
funcbion bo be 

^ Cr> « sin mr 
w ^ 
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In this basis the surface potential from equation C2.3.13> 
can be written as: 


KCr ,r' ,U . ,E> 

B O coh 


»er 




h > 
o » 

h Cxi" > 
o e 


Y Cr> Y CrO 

oo oo 






il/2 


where, x ■ 2C U<:E>-E > j and L ■ Cl,m> ■ <:0,0>. 


Substituting the Hankel function, = -Cl/ix> expC-ix> and 

Bessel function J Cx> «» sin x ,where x = xr, we obtain 

O — - 

X 


KCr ,E> 

& s 


1-1 «r 


2r 


Y cr> Y cr > 
oo oo 


C3.1.3> 


Now, in the ^iven basis and with this surf ace potential, the 
effective Hamiltonian matrix element is ^iven by 



+ sin nr sin mr 

d 8 



8 


C3.1.4> 


where the overlap matrix O is riven as 

nm 
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s;ln<n-jn>r' sln<n4-m>i> 

fi £ 


Cn-jn> 


<n4*m> 


<3.1.5a> 

Fojp a slngl© TnaiTii>©x* bassist Jfuncbion, wban in « bba above 

mabrix element will become 


H 


C J + V_Cr» O 

^ 1 r»n 


+ — sin 2nj> 


and 


sin nr 


i( 


^ , r - sin 2nr 
2 I s a 


2n 


<3.1.5b> 


The unknown x in H can be obtained from CPA equation given by 


c Q <V ,U,E> + a-c> G <V ,U,E> * G CE,U>. 

nr. A nn B r.n 




Equivalently this can also be reduced to 


I IB O I 


r A o vo^ '1 

1 nn nn 

+ 

nn nn“ nn | 

1 

2 

2 

j 


J 


- - IB O 

nn r 


CE+V> H 


+ cv V +EV > -A CV+E> 

A B nn nn 


C3.1.7a> 
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where we have used the following relat^ions : 


V « 


V « 
A 

r»r> 


c V Cl“c> V 


A 

B 

a-c> V + 

c V 

A 

B 

« ^ o +5 

sin 2nr 

2 rin 4 

s 

B - 1 

Ceos 2nr 

nn 4 

2 

U ^CE> - 

>E> 

coh 

2 


<3.1.7b> 


Equation<!3.1.7a> is a cubic order polynomial equation in 
We caiculabe bhe densiby of sbates for ail -these complex roo-ts of 
the polynomial. The negative density of states with Her^lotz: 
analytic properties is the corx’ect choice of the density of 
states. 


3.1,1: RESULTS 

Fii^ure 3.1 Ca> shows the density of states for » -1.0 au 

and V ts -5.0 au .r « 2 au >c « 0.5. The broken curve is 

fi B 

obtained for n « 1 and the full curve for n » S.The essential 
qualitative f eatures are already reproduced for n ■= 1> whereas 

increasing n leads to shifts in the peaks as well as in the 
broadening. Essential differences occur at higher energies. 

Figure 3.1Cb> shows the result for the same potentials but for 
concentration c « 0.9 at higher energy regimes for n “ 4 Cbroken 
curve> and n « 8 Cfull curve>. The two differ very little except 
at very high energies and even then the ratio of the difference 
to the average peak height is only about 4%. 




Figure 3.1 : Ca? Densities of states for V^— 1 . 0 


au,V =-5.0au,r =2 au and c-0. 5: .single member 

basis! eight" member basis. CW Difference in the 

densities of states for a four member 

basisC :> calculation and an eight member basis 

calculationC ^for the same parameters as above 

with c=0. 9 at higher energies. 



Figure 3.2 ; Ctensities of states for s states in 

fiat potential binary alloy with CsO V^=-0. 5 au and 

V =-l . 0 au.Cb? V =-1.0 au and V =-5.0 au and CO 

Bab 

V =-5. 0 au and V^=»0. 5 au. For CaDand 

CbD:C )c=0. 3;(-*-)c=0. 5;( )c=0. 7; for Cc^-A lc=0. 1 ; 
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To lllus’tr'at^e "the effect of' t»he concent.X'a'tion and ly ~V 1 

* A B * 

Cthe two parameters characterising disordered alloys>, we show in 
figure 3.2 three sets of' alloys tor n * 1 ^broken curves^ and n ■* 
8 CfuH ciirves> and c * 0.3, 0.5, 0.7. 

Figure 3.2Ca> is for V * -1.0 au and V * -1.0 au. This is 

B A 

the overlapping-band weak disorder case. The peaks corresponding 
to the constituents overlap to form a broad structure; the 
relative weights of each constituent are reflected in the 
concentration. 

Figure 3.2Cb> is for V = -5.0 au and V == -1.0 au. This is 

B A 

the s^plit-hond strong-scattering regime. Each constituent 

contributes distinguishable structures whose weights again are 
reflected in the concentrations. 

Figure 3.2Cc> is for V * -5.0 au and V » -0.5 au and c = 

B A 

0.1 and 0.9. This is the so called impur'ity-banct 

strong- scat-tering regime. It is here that, from otir experience 
with tight-binding calculations, cluster effects Cnot reflected 
in the CPA calculation here> are expected to dominate. 

3.2: ELECTRON WITH D-SYMMETRY IN SQUARE WELL POTENTIAL 

Now we shall look at the density of states for d-like states 
in a flat potential. Our basis function is then: 
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■[( 


CrO • I I 3/nr,i* - 1/wr I sin mr* 


J sin mr* - cos wi'j 


cos wi' I Y <n> 
22 


C3.2.1> 


The surface potential in this case will be 


KCr ,U,E> = C 1/2 r^> 

B B * 


-i >e^r^-4ae^a>^-*91>£r +P 

S 6 S 

2 2 

r -*31^er +3 




Y Cr> Y Cr’> 
22 22 


<3.2.2> 


where we have taken the following Hankel function 


. ^ . A f ^ ^1 

h Cx> *i - + — - — e 

2 X 2 3 

»- X X 


and the Bessel function Is 


J Cx> ■= r - - - 1 sin X - — 
L X J 


COS X 


with X s= xr. 


Here the effective Hamiltonian matrix for a single member 
basis with m = n is 


H ■ Cl/2 n^ + V^<:r> > O + C + B 

nn I nri nn nri 


<3.2.3> 


where the overlap matrix is 


nn 


1 

2 n 


r 

2 I B 


sinanr 1 3 cos2nr 3 cos2nr 

- a l-t- r s “"3^ s 


2n 


3 

n r 

fi 


n r 


3sin2nr 


n r 2n r 2n r 

s s fi 


] 


C3.2.4> 
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nn 


[ f - — 

+ 3 

+ 1 — 1 

L 1 

4 3 

n r 

2 J 

nr 


sln^ra* 


r 54 

6 

+ ^ 1 

{^3 4 

'"nr 

B 

3 2 

n r 

s 

J 

-L 

9 

3 1 

1 

4 3 

n r 

B 

nr 

fi 


cos nr sin n 2 > 


cos nr 


B 


nn 


Mr * j_ 1 

2|l tf!5 43 2 J 

nr nr nr-^ 


sin^nr 


nr nr nr 

St B B 


- f ^ ,Jl-] 

I 5 4 3 2 J 


COS nr sin nr 

3 2 1 S B 

nr nr 

S B 


[ 4-3 ] ] 

n r J 


and 




B & B 

-“«^r^-“3i«r +3 

S B 


Hero CPA equal. ion reduces t.o t.lie jfoJlowing St^h order polynomial 
equatrion in ae. 


iB r 

nri s 5 
n 


2 2 
A r VO r ^^2 

nn a nn a - tJ 

— — B , riT’.- 


] ”* 


3ir A 

a nn 


3ir V O 

a nn 


^^^a®nn - ir® B E -iV B r® 1 x® 

— nn B J 


B r»r» 
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-( 


EA r + nn - nn -VO r E - “nn - 4B r^ E + A Vr* 

nn » ^ 2 r»r» « — r — ■ nn s nn & 


-V V O - 4i'^B V 

A B nn B B nn 


]«“ 


-r 


31EA r + Sir E V O - 91r- B E + 31A Vr -31r V V O 

nn s B nn s nn nn s s B A nn 


-91r« VB 1 « 

» nn J 


3EA 


3VO E -9B E -3V V O 

nn nn A B nn 


+3A 


V -9B V 

nn 


O 


<3.2.5> 

( n^ *1 ->« — 

- E O + G . V and V are as ^iven in equat^ion 
2 I m“i nn 

C3.1.7b>.The correct, choice of density of st.at.es can be obtained 
according bo bhe criterion given in section 3.1. 


3.2.1: RESULTS 


Figures 3.3 Ca>-Cc> shows the d-state density of states for 
the corresponding parameters of the s-states shown in the earlier 
figures. The structures are sharper. This is to be expected 
since the d— states fall away more sharply from the well centres 
than do the s-states. The qualitative features are well 

reproduced by the single member basis. 




Figure 3.3 ; Densities of states for the d states in 
a flat potential binary alloy with the same 
parameters as the corresponding parts in figure 3. 2 
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3.3: ELECTRONS WITH S-SYMMETRY IN HYDROGEN LIKE POTENTIAL 

To sLudy effect, of Lhe well shape> we next, carry out. a 
similar calculation for a 1/r type of potential. Ve shall look 
at the s like states in the 1/r type of potential :VjCr]>=*CVj>/r 
where V «» V or V . The basis function is 

A B 

3/2 

X Cr> » 2ryi expC-jnr^ Y Cr>, 
rr, 00 

Though t^he shape of -the pofent^ial is differenf but. if sfill 
refains fhe s symmefry in ifs lowest energy state. Here the 
surface potential will be same as given by expression C3.1.P>. 
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IB O 

nn nr* 

2 


H 



V c o * 


nn rin 


- O A 
2 nn nn 


] 


2 




IB G <:e+v> h 

nn nn 


and 


+ CE C A - VE + VC 


r.r. r., - A - V V C > ■ 0 

rtn nn nn A B nn 


C3.3.2> 


A = E n^x- - 2 n* 3 exp<-2nr > 

“ fi * 


Haro also wo Itavo usod t^lio so'L of" rolat^-ions givon in 
eqtiat/ion C3,1.7b!>. This is ag-ain a cubic order polynomial in x 
and bhe choice of* bhe root is “the same as before. Once we geb 
the root^the calculation of local density of state follows the 
saune procedure as described in previous cases. 


3.3.1: RESULTS 

Figure 3.4 shows the effect of the well shape. Ve have 
taken the two potentials to be Vj<r> = CV^^/r where V can be 
either or V^. The qualitative features remain unchanged from 
the flat potentials>although the quantitative change in the 
location and the widths of the structures are apparent. The 
structures are sharper than the flat potentials^ which is a 
reflection of the fact that the flat well s states fall off as 
1/r, whereas the 1/r potential wavef unction decays exponentially. 
The density of states also correspondingly decay at higher 
energies, and the peak widths are narrower. 





ENERGY 


CH) 


Figure 3.4 : Densities of states for the s states in 
a binary alloy with potentials of the form 
VCrD=V/r, with V taking the values and with 
probabilities c and Cl-cD.The parameters are the 
same as in the corresponding parts in figure 3.3 
and 3.3. 
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3.4: ELECTRONS IN NON-SPHERICAL POTENTIAL 

In Lhe previous discussion, we hiave examined several model 
cases embeddinj^ spherically symmetric potential in Coherent 
Jellium. In this section we study single sites and we extend 

the idea of the Embedding method to non—spherlcai potentials and 
study the effect of non-sphericity in the Coherent Jellium model. 
Here we consider the potential in subspace I to be of the form 


VjCr> = C-Vj/r + X r cos^ > 


C3.4.1> 


where the equipotential surfaces are ellipses, figureC3.6>. Here 
X measures the deviation from sphericity, Vj can take the values 
or with probabilities c and <l-c> respectively for a binary 

distribution of the alloy. For small values of X the chosen 
basis will be CLuban & Nudbr-Blum 1977) 


;t<r’) ■ 2. exp<-r) X Cr+1/2 r^) 

<3.4.2) 

The effective Hamiltonian matrix in this single member 
basis, i.e., with n * 1 is given by 


H » 

nn 


A - 

nn 


C -i 

nn 




B 

nn 


C3,4,3) 


where 


= 2n 


- 1/3 


[ 1 - 4^3 


I -2/3 

4 


I -1/2 

5 


I +r 
2 s 


expC-2r > 
s 
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CENTR 




B 


c - 

nr^ 


* 2n expC~2r‘ > 

s 

[4^ X'|l -ij 


'Siser - + 3> 

S BO 


1 - i» r 


and 


expC“2p n/2 I ,n > 2 

c a Fi— 1 


C3.4.4> 


and 


I « -1/2 r expC-2r )-l/4 6xp<:-2r >+1/4 

4 s s » 


The overlap matrix is given by 


-4^3 [l, 41^ *1/4 * I, ] 


The surf ace potential will be 


KCr ,r%U,E> 
a a 


l-i«r j 


Y <r> Y^^CrO 
oo oo 


2r 


1-21 xr - 


1-1 xr 


2 2"^ 
X r j 

— J 


y Cr> Y CrO. 

40 40 


C3.4.5> 


Here we have 


used the foUowlng Bessel and H^el functions 


L'JRARt 

11255 ? 


J . £S^ , expOxl 
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sin X cos X . , ^ 

J/"’- i — , V'" 




expCix> 


and X can b>® calcula'Led from t>he CPA oquat.-ion whdch neducas t-o 
■the foUo'wing fourth order polynomial equation in x. 


C-iB C r^^x** + i<3B C r^-ir EO C +ir C A -ir V >x^ 

nn nn m 2 nn nn c « nn nn & nn nn » nn 


i^CEO C -C A +C V+4iB G Er® + 4iB C Vr^>x^ 
2 nn nn nn nn nn nn nn s nn nn a 


+CiE^O C r -lEC A r +1EC^ Vr -3B G Er^-lEO C Vr +1C A Vr 

nn nTi a nn nn & nn s nn nn s T>n nn a nn nn a 


-iV V C^ r +3B C Vr^>x 

A B nn s nn nn a 


+<:2iBCVr -E^OG-ECA-G^TO-iBCEr -EOGV+GAV-V V C^> - 0 

s 8 A B 


C3.4.6> 


3-4.1: RESULTS 

f 


In this calculation or V with probability c and 1-c 

1 A » 

in the x*eg:ion I with X « .2, 3 and .8 nespectivly. 


3.5Ca!^ shows donslty oT states 
“l.au ,c«,5 and with the above given values 
ovenlapping- band weak* disorder case. The 
to the constituents ovenlap to form a broad 


for V * — 3au>V * 

A B 

of X. This is the 
peaks corresponding 
structure. 


As the 
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deviation fi'Oin sphericity increases overlap increases and the 
peaks are shifted to the lower energy side. 

Figure 3.5 Cb> shows a result for V « 5. au> V * 1. au and 

A B 

c » .5 with the above values of X. This is the split-band strong 
band^ strong scattering case. Here the peaks become wider and 
are shifted towards the low energy. This effect is particularly 
prominent for low energy peaks. 

Figure 3.5Cc) shows a result for V « 5. au > V » 1. au and 
c « .1 where the values are same as taken above. This is the 
impurity band strong scattering region. Here the more prominent 
effect of non-sphericity shows up in the lower energy impurity 
band. 


These results are to be compared to the spherical 
potentials case for the parameter shown in figure 3.4. 


Figure 3.6 shows non-spherical equipotential surfaces. 




Figure 3.5 sCa!) Densities of states for V =-0.5 

A 

au,V^=-1.0 au,c=0.5,ij =a.au ; for 'x = 0.5;-»- 

for X=. 8,-. for X= 0.2. 



ENERGY V n 


Figure 3.5 ;Cb> Densities of ststes for V =50 

A 

au , c=0. 5 with the same X’S as given in 
figure 3. ICaD . 



DENSITY OF STATES 



Figure 3.5 iCci> Densities of states for V =5.0 au> 

A 

=1,0 au c=0. 1 and with X’S as given in figure 
3. ICaD. 





CHAPTER IV 


AN EMBEDDING- CPA PORMAUSM IN ONE DIMENSION 

The KXR method proposed four decades a^o by Korrin^a <1947> 
and Kohn and Rostoker C1954> remains a powerful technique for 
accurate calculation of electronic band structure for pure 
metals. Later this method was adopted by Qyorffy and Stock to 
calculate the electronic band structure of random metallic 
alloys. Butler C1976> developed the one dimensional analog of 
KKR approach. In this chapter we review the one dimensional 
analog of KKR approach and then develop the Embedding-KKR-CPA 
formalism to calculate the electronic density of states for one 
dimensional random chain of flat potential wells. 


4.1: KKR FORMALISM IN ONE DIMENSION 


We consider a one dimensional system with Hamiltonian 


H 


S 2 + UCx> 
dx 


where 


U<x> ■ Zv <x-na> 
n 
n 

C4.1.1> 


and the summation is over n atomic sites with lattice parameter 

a. We assume v Cx— na) to be symmetric and vanishing for r* | x | 

r» 

greater than some muffin-tin radius r^^. We also assume that the 
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dlst,anco between the centres of the potentials is enough, 

so that they do not overlap. If the potential is situated at the 
origin then we have 

V Cx> ■ 0 , r ■ |x| > I- , where r < \/ 7 . a 

mi 


For a single 
C4.1.1> we have 


symmetric potential satisfying equation 


[- 


dx 


+ UCx> 




V/<x> 


C4.1.2> 


We shall expand the wave function about the centre of the 
potential in terms of symmetric and antisymmetric functions, 
analogous to the spherical harmonics used in KKR theory; 


Y Cx> « 1/^2 

Y^Cx> « CI/T 2 > signCx> 

C4.1.3a5 

We shall also need the one dimensional analogs of the 
spherical Bessel functions : 


j^Cz> ■ COSCZ - lfT/2 >, 1 “ 0,1 
n^Cz> ** sinCz - 1» t /2 >, 1«0,1 

<4.1.3b> 

Now the wave function satisfying equation ^4.1.2) can be 
written as 

yv 

V/CX> - 2 E a^ Rj^acr> y^<x> 


<4.1.3c> 
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where G^Cxpt' > Is the free electron Green fxjnction and 


m t"' s 6 + r r t”* G' CR -R 


i 2 1 1 a pi^w i ta 


a a 


C4.15b> 

is the scattering path operator introduced by Gyorffy and Stott. 


Fourier “transiform^ 


any 

arrangement. 

of 

If 

we rest.rict 

our 

a 

lat.t.ice^ t.hen 

t,he 

by 

means of a l3t.t.ice 


^Ck> ■ E E exp<:-ikCm-n>a> ^ 
12 rt\ n 12 


■ \ I * E \ O' ^ac> ^ck> 

112 1 113 92 

a 

C4.1.5c> 

where 

G' Ck>» expC-ikpa> G' Cpa) 

.1, ■’’‘O 

C4.1Sd> 

equation C4.1-5c> is a two by two matrix equation for t which has 
"the solution 

Ct~*> = t'* 6 - G' ,Ck> 


C4.1Se> 


The poles of scattering path operator are determined by the 
determinant equation : detCr > "0. which are also the poles 
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the Green function and thus the ei|^en values of the Hamiltonian. 

Then 

det Ct~* 6^^ - G'^ CIc>]« 0 

C4.1.6> 

Is the one dimensional analog of KKR equation. We are mainly 
interested in the Green function G<xpt' > when x and x' are both 
within the same Wigner-Seitz cell. Let us suppose that x and x' 
are both within the Wigner-Seitz cell at the origin <but in the 
zero potential muffin- tin plateau region >. Then using the 
spherical harmonics expansion of all terms about the origin the 
Green function can be reduced to 


G<:xp<'> - G^Cxpf'> + E ak>"^h^<kr>h^Ckr'>t° Y^ci>Yj^<:xO 


+ E [j Ckr>+<:ik>"*tj^ h^Ckr>lYj^Cx>G®°j^ 
1 2 


J Ckr' >+Cik>"‘t^ hj^<kr' >yj^Cx' >1 
- ^2 2 2 2 ■* 

C4.1.7a> 


where we have used t^he relations 


l l 
1 2 


•■r I ^ 

1 12 


0 -oo , o 

1 I ^ 1 

1 12 2 


and 


E 0' 

I 

2 


oo 

,ac >T - 

I m II 

12 2 3 


G t 
>• ^ 

12 2 


C4.1.7b> 


_ OO 

E I ® I <k„> 
I 1 2 2 9 n 


t.°e 


r. 

i 3 
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and 


-i 


^oo . —1 ^ oo 

^ Ct - t <5 > t 

4 2 i » 2 1 12 '■a 


<4.1.7c> 


For X and x' bot.h at one of the boundaries of the cell at 
the oriel'™* expression <4.1.7a> reduces to 


Gfxpt'> ■ £ e <;r,r'> Y <x>Y Cx' > 

, t Li 


where 


C4.1.8a> 


G Cr,r' > - R. Ckr > 

l I < 

k sin6^ 


tCkr >+ R Ckr > 

I > I > oo 

T 

k sind. 


<4.1.8b> 

We have to calculate , i.e., the scatterine path 

i 2 

operator for a path which leaves from and returns to the orlein. 
From equation C4.1.5c> and equation C4.1.5e> this is 


m CaX2fi>J dk <t"‘-e' Ck»"‘^ 

12 4 ; 


<;4.1.9> 

In one dimension taking only 1-0 and 1 partial waves in 

OO 

the expression the elements of this path operator matrix, 


can be obtained as: follows:: 




OO 

T 

OO 


tan6^ r cos8-f^ ^ 

1-t-tand tan<!> J cos8-x 
o ^ _ 
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• -k i^BnS f cos^f 

^ o da 

l+tran6 J cosa-x 

o i 


where 


00 

T 

01 


•T 


.OO 

iO 


- k 

2n 


+n 


tand t.an6 f 1 sina 

— ^ 1 cosa-x 

1 4''Lan<5 ‘Laii6 J 

o i 


X * cosC^<5^4-6^> 

cosC6 -<S > 
o i 

jT * cos^"tan6 sin^ 
o o 

jf * cos^t?~t»an6^sin^ 
p « ka 


C4.1.10> 

The off-dlaeonel elements - t°° - O, because of the odd 

parity of the integrand. So the scatterinc path operator matrix 

in one dimension is 


Making change 
C4.1.10> and can be 
in “the folio wing way: 


oo 

i 

in the variable z » expa0>, the integrals 
performed integrating around an unit circle 


^OO 

oo 

O 


^oo 

11 


+ tT 

1 r cosa-f ^ 
2fir J cosa-x 


1 r ^ z^"2fas’H 
2ni J z z*-2xz+l 


-n 





6 


Tlie of" t/his ssimple poles > 2»0> 2*2 end 

2**2 whene 2 2 "I. For* x noel end less t*hen 1 bot^li 2 end 2 

7- ^ > < > 

will be on the unit clncle. However if we let E have a positive 

imaginary part^ one of" the pair between z and z will move 

< > 

outside and the other one inside the unit circle to satisf"y the 
relation 1. Thus equation becomes 


I 


1 ± 


Cx-^O 
Cx ~1> 


The si^n is uniquely determined by ^ivin^ E a small positive 
imaginary part and seeing which of the two roots 2 ^ or 2 ^ lies 
inside the unit circle in complex energy plane. The path 
operators will become 


T 


00 

00 


-k 


tan6 
o 

1-K,an6 tan6 
o 1 


ICx^f > 
1 


T 


00 

il 


-k 


tan6 

i 


l+tan<5 tan6 
o ± 


ICx,f > 
o 


<4-1.12> 


4.2: KKR-CPA FORMAUSM IN ONE DIMENSION 

Thus the one dimensional analo^^ of KXR formalism developed 
by Bubler gives the Green function for a single scatterer. When 
we apply Embedding method to one dimensional arrays of flat 
potentials subspace I wiU be a periodic system of identical 
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sca-tt.«x*ers. I-t was sueeest.ed by Anderson and McMillan C1967>, 
for liquid mebals, that the Green function within the 
Wlu'ner'Seitz cell at the orlg’in for a periodic system is equaJ. to 
that for a single scatterer in a uniform medium. 

According to the Embedding theory, subspace I in one 
dimensional arrays of flat potentials model can be tx^ated as a 
single impurity problem in an effective medium and the Green 
function with Neumann boundary condition at the interface S of 
the two subspaces I and II can be easily obtained. 


The CPA can be derived requiring that the Green function for 

a system having scattering amplitudes t and t at the origin and 

A S 

iuhe coherent^ scattering amplitude t^ at all other sites be equal 
on the average to the Green function for a system having t^ on 
all sites. Equivalently this can be written as 


c 

A 




oo 

C,A 


+ C 

B 


oo 

T 

C,B 


T 


oo 

C 


where is given by equation C4.1.12> with tan<5^ 

c 


C4.2.1> 

tan6^ and 


C , A,li C , it A,i C,i 

Substituting one into the other yields an equation for 

on- energy shell t matrix, t : 

c 



c t 

A 


B B,l c,l A, I c,ll c,l 
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The inverse of the on-ener^y shell t matrix for a one 
dimensional array of flat potentials can be calculated to be of 
th& form 


-k tcot6 - a. 


C4J2.3> 


Eqxxat^ion ^4-2.2> is an adl^ernafive form of "tHe CPA aquafion. 
for random binary aUoys. To solve this self-consistently a 
suiteJble initial choice of t = c t c t is made and 

C,l A AA B B,l 

usinc in equation C4.2.3> we can calculate the phase shifts and 

substituting in equation C4.1.12> an initial value of can 

c ,ll 

be calculated. Substituting these values in equation <4-2.2> we 

obtain a new value of t ^ . This goes iteratively until 
converges. Now we can calculate our Green function of subspace I 
at the interface, substituting the above Rj^Ckr> and j^Ckr> 

are defined in equation C4.1.3d> and <4.1.3b> whereas the phase 
dependent terms are obtained from 6^*^. 

4.3 ! EMBEDDING - KKR- CPA FORMALISM IN ONE DIMENSION 

Here we are going to develop an Embedding-KKR-CPA theory in 
one dimension. The one dimensional space will be divided into 
two subspaces, i.e., I and II with an interface S at |x|*«r^ as 
shown in flgureC4.1>. 

Where *Cx> is the eigen function of the Hamiltonian H of the 
whole one dimensional space satisfying the Schrodinger equation. 
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H »Cx> ■ E »Cx> 


wi'th energy eigenvalue E. 


C4.3.1> 


To solve t-he Schriidlnger Equa-tion in -Lhe region I we find 
■the expectation value of the energy of the wave function ^x> of 
the full space, defined in svibspace 1 as a trial wave function 
tpCic'y and in the subspace II as the exact solvition v'^x> of the 
Schrodinger equation at same energy E with an extra condition 
that the solution matches correctly on the boundary surface S : 


and 


2 

^2 + U<x>-Ejv'<x> 


C4.3^a> 


V/Cx > «* i^Cx > 

d S 

<4.3^b> 


The ener-^y expecbablon value of “the Hamil-tonian of t-he full 
space is 

CO 

J dx ’|f*Cx>H'Kx> 

^ 00 

E * 

00 

J* dx 4^* <x>'®<x> 

“00 

<4.3.3) 


Because of the discontinuity of the derivative of the wave 
function at the Interface S we take an c neighbourhood of the 
interface S to evaluate the integrals in subspace I of equation 
C4.3.3) and then take the limit e to be zero as shown in the 


figure C4.2>. 
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Figure 4.1 : Partition of one dimensions] 
into two subspaces, t.e, I & II. 
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Figure 4.2 : Epsilon neighbourhood of subspace 
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Thus subst.lt,ut.ln« eqti£it,ion C4.3.2a> in e<iuatlon C4.3.3) we 
obt^ain : 


-r ~£ r *£ 

E J dx «« Cx>4'Cx> + E J dx 4^ Cx>4f<:x> +E J dx 4'*Cx>4'Cx> 

“00 “r -£ I. +£ 

m fi 


-l' • 

6 



-00 


dx 4'^Cx>4'Cx>+ 


-P 4-^r X* — £> 

* * * * 

J dx «< Cx>H«<;x> + J dx ^ <:x>H<^Cx> 

-P -P 4-r 


P 

« ® * 

+ J dx ^ Cx>H«^Cx>+ E 

p 

a 


J 


00 


dx ^ Cx>^<]x> 


p 

3 


which can be reduced bo 

p -►£: -p p 

® a * a * 

E J dx 4^ <x>4'Cx> - J dx 4* Cx>H4'Cx> + J dx ^ <x>H^<:x> 

3 ~P •P 

fi a 


P 4-^ 

* 

+ J dx Cx>H»Cx> 
p -e 


C4.3.4> 


Leb us now evaluabe each inbe^pal sepapabely. The f'ipsb berm 
wlU be 


-p +s -p ^ 

fi fir* 

J dx 4f*Cx>H4'Cx> - J dx 4-* Cx>H^j4<x> 


~p -e 

a 


P -4? 
S 
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-r +£ 

-»*Cx> |^'l'Cx>| “ +J dx Cx> |^>PCx> 

~P 

« -P “iP 


With limit £ going to 0 , equation 4.3.3 becomes 


" 1 -I. -/<-i'^>^V'Cx> I 1 

>- * s a-* 


Similarly t^he -third term will be 


r 

s 


J dx «<*Cx> J-^2j4'Cx> - - i3P*Cx>|^4<x> 


r 

a 


r -ts 


r 

s 


r 

* j ^ j 

+ r dx 5 ^ <:x>2 ipcx> 

J dx dx 


r -s 
s 


-t 


A A J 


Now in the limit £ going to 0, substituting these two in 
equation C4.3.4> we obtain 


Ej dx 4 >* C.x:><pCx:> m J dx <^*<x>H<^Cx> 


-r 


-r 


L fi »-* 


-p 




C4.3.5> 
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According to Embedding theory now we express the derivative 

In subspace II In terms of the trial wave function 

■ * II 

<^Cr > in subspace I by using the Green function G Cxpc' > for 

region II » satisfying the equation 




[“ix* + U<x>-Ej6“cxpt'> - «5<x-x'> 

C4.3.6> 

where x and x' are in region II, i.e., |x| >r^. Multiplying 
equation C4.3.2a> by G^^Cxpt' > and equation C4.3.d> by and 

substituting one from the other we get 


"dx^ \Kx>1g^^Cx,x' > + 



G^^Cxpt' > 


lyCx> 


-6Cx-x' >v'^x> 


<4.3.7) 

Integrating through the region jx] > r^, exchanging x and x' 
then substituting the boundary condition on the surface the above 

equation will reduces to 


-r -e 

V-Cx>- G^^Cx'pc 

^ -00 


+ G^^Cx' ,x>^,V'Cx' > 


oo 




With limit £ going to zero this equation becomes 




^ ^4 4 n.+ <.T.face of the two subspaces will 

So the wave function at interlace oi 


be 
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V^Cr > - <^<:r > ■ G^^C-r ,r >^^,V'Cx'>| - G^^Ci* ^ .y/Cx' > 

^ a a a a dn'’^ l-r a a dn'^ 


a a 

C4-3.8a:> 


V/C 


-r >■ (^<-r >■ G^^C-r ,-r* .y/<x'>| - G^^Ci* ,-r V'Oc'il 

s s SB dn ' -I* a a dn 


C4.3.8b> 


JI. 


Mulblplylnc equation <[4.3.8a> by 6 C-r* ,-i* > and equation 

a a 

<;4.3.8b> by G^^C-r ,r* > and stittnactin^ one fx-om the othei* we 

a a 


ob'balii 


(*<r >G^^C-r ,-r >-<^C-r- MS^^C-r ,j> > 

Q a a e aa 


- ^G”<r^,-r.^>G^^C-r^,r^>-a^^Cr^,r^>G^C-i-^,-i-^>]|^,V'<x' >1^^ 

C4.3.9> 

A^ain substituting D for the bracketed term on the right 
hand side of equation C4.3.9> we get 




3«C-r ,r >] 

8 S J 

<4.3.10a> 

Similarly multiplying equation C4.3.8a> by G <:r^,-r^> and 
equation C4.3.8b> by G^^^^^s'^s^ and substituting one into the 

other we get 






C4.3.10b> 

substituting these two expressions in equationC4.3S> we get 
Ej *dx - J dx Cx>H^Cx> 


-r 


-r 


L 0 a * 
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p^Cr.^ -r^>^Cr^>-Q^^<r^,r.^><?i<-i.^> j j 

C 4 . 3 . 11 > 

This expression chives us E purely in terms of a trial 
function <^Cx> defined in subspace I and on the interface S. The 

variational minimization of E with respect to the trial wave 

e 

function 4> Cx!> t.e., 

‘SE 

__ . 0 
Sip Cx> 


^ives 


r 

E 0<x> » 155^^ <^^x> 


— < 


6Cr»'Hi^> ^ <;6Cx> - 6Cn“n^> ^ ^Cx>l 

a a J 


6Cn-n > J fo^^C-r ,-r >-G^^C-r ,r >«^C-r >1 

aDL aa^a ea aj 

- (5Cn-Hn > r fo^^Cr ,-r >^<r >-G^^Cr ,r >i^C-r > 1 

aDL a'a^a aa aJ 


C 4 . 3 . 12 > 

This is the one dimensional analog of Embedded SchriSdlnger 
equation that acts on the subspace I alone and the effect of 
stit>s:p 3 C 4 B II appe 3 r*s 3 S t»h© sui'f' 3 C<B pot/ent^ial opei* 3 t»oi> 


j6Cn-n^> i GC-r^,-r^>-6Ci»+«^> 5 


G<-r ,r >+< 5 Ci»+« > s I 

s' a a D * * J 
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THo t#©r’Tns2 <5(nHhn ) ^nd <5Cn—rx ^Cd/dn > enstut^^ trhat^ 'LHe 

mm mm 

effect^lvo Hamllt-onlan remains Hermitian In region I alone. 

4.3.1: RESOLVENT OPERATOR IN ONE DIMENSION 


Fop acLiial calcula'tion of one dimensional pesolvenf opepa'top 
we need a peppesenfafion of effective Hamil'tonian in a stiiLable 
bausis. So we expand oup -trial wave function in subspace I in 
terms of a countable basis ;^Cx> which spans the subspace 1 and 
reflects the symmetry of the solution we wish to study in I. 


0Cx>=E 

i 


C4.3.13> 

Usinjg^ this basis in ecfuation C^.S.ISD and imnimising^ E with 

respect to variation in the effective Schrodin^ep equation in 

one-dimension can be written as 

E h a *= E r O a 
Ik k f- Ik k 

k k 

C4.3.14> 


where 


p 2 


-P 


-**C-r^>|^ [o^ ' Ci-_ ,-1-^ 


C4.3.14> 
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and "th© overlap ma-trix is 

r 

* ^ 

^ik“ J 

* r 

fi 

C4.3.15> 

Here t^he Green fundion in subspace II can be expanded as 


G^^Cxpc' >-2 G^^Cr,r' >Y^<x>Y^Cx' > 


C4.3.16a> 


where 


Gj^Cr,r' > 


1 

B 

ksin6. 




<kr^>- 


R'Ckr > „ ., 

I a P^Ckrj^> J 

F'Ckr > 

I a 


C4.3.16b> 


and 

F^Ckr>* |^jj^<kr>+ j 

k sin6j^ 

C4.3.16c> 

Now analogous to the three dimensional embedded theory the 
matrix representation of the embedded Green function or resolvent 
operator in one dimension will be 


G 


Ik 


<H-EO>‘f; 

Ik 


C4.3.17> 


where the local density of state is defined in ch^ter II 


CHAPTER V 


ONE DIMENSIONAL MODELS 

In t-hls chiapLer' we shell illust/Z'et>e otur* one dimensionel 
analog of Embedding formalism taking single and double muffin-tin 
potential wells of the binary alloy type with probabilities c and 
<l-c> and potentials Vj<r> - -V^ or -V^ in the region I. In 

region II we replace translationally symmetric muffin— tin 
potential of the binary alloy type by an effective CPA medium and 
evaluate the path operator and siirface Creen function using 
KKR-GPA theory, Butler C1976>. Let us now calculate the locad 
electron density of states embedding single and pair muffin- tin 
potential in the one dimensional effective CPA medium. 


5.1 : SINGLE MUFFIN-TIN POTENTIAL IN SUBSPACE I 

Ve choose our single member basis function to be of the form 

XCxy « 2C1/N> Ca ;^:Cx)+a^;t<^x» in subspace I, where N is the 

i<5 , i.S 

normalization factor; a and a are given by e o and e i 

O X 

respectively f or 1*0 and 1 partial wave states. Expanding 
;*:^Cx) and explicitly in terms of equations C4.1.3a>, 

C4.1.3b>, C4.1.3d> we obtain 


with 


X <.x> ■ R Cr> Y Cx) 
^o o o 


R <r> * cosCkr+6 > 
o o 


r < r < r 
o s 

-r < r < r 
o o 


■ cos xr 
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and 


■ cos Cki*-6 > 
o 

y Cx> » R Cr> Y Cx> 
i 11 


-I* < r < -r 

a o 


with 


R Ci'> » sin Ckr'+6 > 
1 1 

» sin xr 

■ -si n xr 

» sin Ckr*-6 > 

1 


r < r < r 

O S 

0 < I* < p 

o 

-- I* < r < 0 

o 

* r < r < 

s o 


C5-1.2> 


with • < E-V Cr», where V-Cr> can either -V or -V and » 

I I A B 

E The phase shifts for 1»0 and 1*1 partial wave states can be 
calculated respectively from the relations given by equation 
C5.1.3> : 


tan 6 Cu tan «r -k tan kr >Ck tan kr tan xr > 

0 o o o o 

tan 6 ■ Ck tan xr - x tan xr > Ck tan xr tan xr + x > ^ 

1 o o o o 

C5.1.3> 

The normalisation factor can be obtained from 


r 

^ J dx ;t:*Cx> ;»:<x> 

-r 

B 


In this basis the overlap element is 


O 




dx ^ » 


-r 

a 


<S.1.4> 


C5.1.5> 


where we have 


calculated these integrals numericaUy using 
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comput.er. Th« Hamilt^onian is 


H ■ E O + <4k/^> I- cos Ckr +6 > sin Cki* +5 > 

L * o so 


+ cos <lcf +6 > sin Oo* +6 > 

Si Si 


] 


+ C4yN*> 1^ cos^kr^-t-6^> ^ sln^<kr^-^5^ > -j 


6 Cr ,r > 

o a 


6 Cr > 
i a a 


C5.1.6> 

where t.he Green function of the effective CPA medium on the 
Interface S is given by equations C4.3.16b> and C4.3.16c> and the 
path operator by the equation C4.1.12> Giving a small positive 

Imaginary part C 10 to the energy we get analytic density of 
states at all energies except at zero and a few other values of 
energies, because our path operator calculation fails to work at 
these points. 


The local density of states is 


n<r.E> 



O 


A 

EO -H 


A A 


C l-c> 

n 


Im 


EO -H 


<5.1.7> 


5.1.1: RESULTS 


In order to Identify the structtxres in the negative energy 
range of the density of states of C10% —90%^ composition of the 
binary alloy, we compare the local density of states of the alloy 
5.1 Ca> with the local density of states in A and B 


as in figure 
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at/Oms in alloy as: shown in fi|fui>e 5.1Cb!> and 5.1<c> 

respocrtlvely. We have taken V ■ -0.3 au» V ■ -0.5 au. Ve 

identify A^, as peaks associated with the A constituents as in 

ficure 5.1 Cb> and with the B constituents as in figure 

5.1Cc>. 


There is some overlap between the two, leading to widening of 

the structures, to form a three peaked band. Peaks A and B are 

1 £ 

wide peaks whereas A is a sharp one. The peak A sits on a 

narrow band, coming from the medium. This is separated from the 
previous structure by a narrow gap. The positive energy bands 
have long widths, so that their smoother structures are not 

identifiable. The characteristic Van Hove singularity structures 
Cwidened by dlsorder> are well describable. 

We analyse the <50% - 50% > composition of the negative 

energy density of states of ICPA, shown in figure 5.2Ca> 
comparing with the density of states of A and B constituents 
shown in figure 5.2Cb> and 5.2Cc> respectively. Ve find there 
are three structures. The left most A^ and right most A^ 

structures are from A constituents of the alloy whereas the 
single peak of B constituent Is B^ situated in between and is 
sharper than A^ and A^. All the negative energy structures are 
well separated. If we analyse the peaks of 'B constituent we find 
there is a sharp delta function Uke peak and a small sharp peak. 
In between these two peaks the density of states does not 
converge because our path operator calculation fail to work at 
these energy points. 










1 -ac 

ENERGY 

Figure 5.1c :1CPA, negative energy partial 
densities of states of the B constituent of the 
binary alloy for V^=-0.5 au. 










DENSITY OF STATES 
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Figure 5.2d:]CPA, positive energy densities of 
stQies for V au«V =-0,5 au and c=0.5. 
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It is pertinent to mention here that the method Tor choosing 
the roots in t he calculation of the path operator 

occasio nall y falls in the negative energy region. 'Whereas f or 
positive energies> one can always take the complex conjugate 
Herjjlotz alternative, this is not possible for negative ener^^es. 
We have no way to obtain the correct path operator at these 
points. Fortunately, these energies are few and widely 

separated. In more realistic applications, the path operator is 

obtained numerically via Ewald's method and such difficulties do 
not arise. Wherever we have such difficulties we have indicated 

this with red arrow. 


On the positive energy side we get as before, bands with 
smoother structures. 

5.2: A PAIR MUFFIN-TIN POTENTIAL IN SUBSPACE I 

Here we replace the region I by a pair of muffin— tin 
potentials of binary alloy types, i.e., VjCr> « -V^ or with 

probabilities c and Cl-c> respectively. We embed this in 

effective CPA medium described in section 5.1. We choose the 
member of the basis functions belonging to each well as 

XM . 2 i *,,0**’ * “i ^ 


2 <:x> - 2 S ‘ + K ’ 

2 


C5.2.1> 
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where x^Cx> is t,he basis of the weU with potential V^Cr> and 
^^Cx> with potential V^Cr>, as shown in fi|;ure 5.2. Here and 


N are the normalization factors and a 
2 o 


expCi6 a 

i 


expCi<5^^^>, b^** expCi6^^^>, b^ ■ i expC16^^>. We assume that 

the total wave function is y/Cxy» ^ * B x • Expanding trial 

wave functions in terms of equations C4.1.3a, 4.1.3b, 4.1.3d> we 

obtain 

V ^Cx> - R Cr> Y Cx> 

4,0 1,0 1,0 

with 


with 


R Cr) » cos <kr-6 > -zr < r < -Cr +r > 

1^0 6 B o 

« cos ae r -Cr +r > < r < -Cr -r > 

A SO BO 

« COS Ckr+6 > -Cr -r > < r < 0 
i^o 8 o 

C5.2.2> 


2>o 


Cx> » R Cr> 

2,0 


Y Cx> 

2,0 


R Cr> * sin Ckr-6 > 

2,0 2,0 


* sin ae r 
2 

» sin Clcr+6 > 

2,0 


0 < r* < Cr -I* > 
8 o 

Ci* -I" > < r < <1^ +1' > 
8 0 so 

Cr* +r > < P < 21^ 
so a 


C5.2.3> 

Similarly we can write the wave functions for 1»1 partial waves. 

Here - CE-V >, - CE-V >, k* » E. The phase shifts of these 

1 12 2 

two wells for 1 ■ 0,1 partial waves can be obtained from the 

r*el 2 it#lons given In equa'tions C5.2.4> : 


,4i 
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* A. - ^ Cr + r > > 

t^ain QWi iso so 

i,o 

< k + 3ie “tan kCi> +r» > tan ae <i* -hr » 

i so iso 


-1 


^ e C k tan ae Cr -hr >-3e tan Ci> + r > > 

tan 6 m isoi so 

i^i 

( ae k tan kCr -hr > tan a <j> +r » 
i so i s o 


^ _ < ae tan a€ <r -hr >-k tan Ci^ + r > > 

tan 6 m 2 2 sO so 

2,0 

C k + ai tan kCr +r > tan h Cr 

2 SO 2 s O 


, C k t-an j« Cr +r >-x tan Cr + r > > 

tan 6“ z e o 2 a o 

2,1 

C K + k tan kCr +r > tan x Cr +r- . » 

2 SO 2 S u 


In this basis the Hamiltonian has the 2x2 
pepnesentation with elements 


H ^EO 'i-'^2 Cx> 

ii ii — r— ** ^ ^ 

PT P -r 
i s o 

- rcc«<2ki- slnCJkl- «°>- oos(:2kr «;> rfnakr^-hsS] 


r cos*C 2 kr^+i 5 ^> ^ sin C 2 kr^+< 5 ^> ^ 

i e Czr , 21 - > 0 ^ 

oa ® IS" 


2 

N* 


H - E 
12 12 


-I* +r 
a o 


''t j • “ dx x,<xyxOi> 


"A ■'■.'-o 


•i 

“1 

CS.2.4> 

matpix 
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2 k r * 

■ N N 1% 
12*- 


b co8C21cr +6 > sinC2kj> +6”> 

‘ O « 1 • 2 


- b a* cosC2ki' ■h5*> sinC2kr- + 6 ^ 


~ irm cosC2kr co8C2la> +6°> 

NNloo si «2 

i 2 

G Czr ,21^ > 
o a m 


+ b a Gos<2kx* 

s : 


G <21* ,2r 
1 s s 


‘> slnC2ki‘ +6*> 1 
2 « *■ 

,2r > J 


V r +r ^ 

H ■» EO ■+• 2 r « o ^ <x>v Cx> 

21 21 j 2 1 

N N r -r 

12 SO 


2 k r , 

N N 

4 2 L 


cosC2kr + 6 ^') sinC21cr 

e 2 a 1 


b*a cosC2kr +(5*i> sln<2kr +6* 
11 8 1 8 2 


— „ a b cosC21cr +<5 > cos<2ki» +<5 > 

NNoo ai 82 

1 2 — 

0 C2r ,2r > 

0 8 8 


+ b*a cos<2kr +6^> slnC21a> +6S 
1 1 a 2 a *■ 

0 C2r > 

18 8 


H ■ E O + 
22 22 


J 


-r +1* 


dx 


N* -r -r 

2 B 0 


- ^ rcosC2kr i-6°> slnC2kr. +5°>- cos<2kr -hSS sinC2kr.^+6Sl 
J^2 L 8 2 *2 * J 
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cos Czkjy -t "6 > ^ sin Cikr* +6 > 


i 


» o 


s i 


G Czr* ,2r* > 

O ft A 


6 Car •ar > 

1 8ft 


and similarly Ihe overlap matrix elements are 


C5.2.5> 


2r 

O - <4/N^ > f * 

Ai. 1 J 


-2r 


dr X 


i 2 


- J ”* dr **Cv> *,Cx> 


-2r 


2r 


o - C4/N N > f . ^ 

2i 2 1 J di' A::2 Cx> a;^cx> 


O - C4/N^> f 
22 2 J 


-2r 


2r 


2 J dr ;t^Cx> 

-2r 

ft 


C5.2.6> 


where t,he normaHsat-ion factors can be calculated from equation 


C5.2.7) 


2r 


J dr 


-2r 


2r 


J dr a:*Cx> 


-2r 


C5.2,7> 


"’'V 
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The calculations of the Integrals are done according to the 
criterion ^iven in section 5.1. 


The Green function of the effective medium on the interface 
S is given by equations C4.3.16b, 4.3.16c > and the path operator 
by equation C4.1.12>. The density of states calculation is done 
according to the procedure given in section 5.1. Since V^<:r> 

and V^Cr> both have binary distributions, we have In all 2^«4 
configurations, i.e., CAA,AB,BA,BB>. Thvis totally there are 4 
sets of the 2x2 Hamiltonian and overlap matrix elements. 
Averaging over the configurations the density of states will be 


n<r,E> ■ c^ Im ^aa 


a-c> 


Im 


O 


BB 


EO -H 

AA AA 


EO -H 

BB BB 


+ cCl-c> Im 


o 

o 

AB 

BA 

EO -H 

EO 

AB AB 

BA BA 


C5.2.8> 

5.2.1 : RESULTS 

The negative density of states structures in 2CPA of 10%-90?i 

compositions are in shown figure 5.3. Here aU the CPA peaks of 

the same constituents splits into an anti-bonding and a bonding 

peak. This is the characteristic of 2CPA <Ahmed & Mookerjee 

1990>. The B peak of ICPA splits into B b and B^a bonding and 

2 

anti-bonding peaks whereas the peak into A^b and A^a 

respectively. The anti-bonding and bonding peaks are not well 
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sep3i*3t>6d. Ttie bonding sbi'vict.ux'e r»f ♦ * 

® or bhese two peaks overlap to 

form > wide structure simlUr to the «,ti-bondlns rtructure, they 

ere end structures respectively. One cen dlstlneuish the 

overlep of the structures in O, .l.e.. bondin. peeks 

of end but In we cen not dlstini:uish between the two 

entl-bondinc structures CA^^ , B^> of A^ end B^. In this cese 

ell structures together form e single bend with humps end peeks 

In between. 


The positive energy structure is quite similar to the ICPA 
case. It does not contain any extra Interesting struct\ire. So 
we have not drawn it separately. 

We analyse the negative energy peaks of C509S-50%> 

composition of 2CPA in figure CS.4a5 in comparison with figures 

C5.4b> and C5.4c> and C5.4d>. Here all the bonding and 

®*^tl-bondlng peaks of AB, BB, AA overlap to form a wide band. 

The sharp peak C comes from the contribution of all these AA, 

z 

ABp BB isharp peaks: which are deno-ted by ^ 2 ^' ^ 2 ^ filjures 

C5.4d, 5,4b, 5.4c> respecblvely. Then t-he st-rucbure drop slowly 

with a small hump G which comes mainly from C d and C b peaks at 

S 9 3 

AA and AB constituents as shown in figures C5.4d, 5.4b> 

respectively. The right most peak comes fiHjm the C^b peak of 

the AB constituent as shown in figure 5.4b whereas the peak 

comes from the contribution of all AA, AB, BB constituents which 

are denoted by G d, G b and G c as shown in figures C5.4d, 5.4b, 

11 1 

5.4c> respectively, forming a wide structure. 
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Here 
st.at.es is 


also t.he change in t.he posit.ive energy density of 
not prominent. 










Figure 5.3 :2CPA, negative energy densities 
states for V =-0.'5 au,V =-0.5 au , and c=0.1. 
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Figure 5.4a jZCPA, negative energy densities of 
states for =-0.5 au.V =-0.5 au, and c = 0.5. 




density of states 



Figure 5.4c :2CPA, negative energy densities o 
states of BB constituent for V^=-0 . 5au. V^=-0 , 3 au 
and c=0 . ? - 




CHAPTER VI 


CONCLUDING REMARKS 

A s:u]>f'^ce pot-ont-lal is dex^iv^d which catn «sidded trO t>hs 
Schr*odin| 7 ep equua'tion Tot 21 limit/ed i^egion of" sp 3 ce|. 1, ‘ho omhed i‘t 
inho 3 subshrahe. This pohenhial^ which is energy dependenh and 
non locals can he f’ound f"3>oin hhe Gr^een f'tinc'hion for hhe huUc 
subs'hrahe. The resulhs are based on a ‘varia'tionad princdple 
which ^ives hhe energ:y of a shahe in herins of hhe wave function 
in I. The embedded Schrodinger ecfuation can be solved by a basis 
set expansion^ for the wave functions of discrete states and the 
Green function in the continuum " . 


— CInglesfield 1981> 

Our work is the alloy generalisation of the embedding ideas 
introduced by Inglesfield C1981>. We have applied this formalism 
to calculate the local electron density of states of Coherent 
Jellium models^ embedding different single muffin^tin potential 
wells in smaller space I, which is the muffin-tin sphere. For 
each type of potential we have done the calculation for three 
different concentrations of the alloy. 

This formalism is developed in the most generalised form to 
take into account clusters of wells as well. Ahmed & Mooker jee 
C1989> applied this methodology to Gohea^ent Jellium models^ 
embedding a pair of miiffln-Lln potential wells. It is 


lnt.eres<,lnc to note that the problem is very similar to the 

problem of a hydro|;en like molecule in the subspace I immersed In 
a Coherent JelUum medium in subspace II. The various integrals 
Involved in this case are the famiUar two centre integrals in 
the finite subspace I. 

The parameters they have chosen are the same as those 

reported by Mookerjee and Bardhan C1989> for one potential 
calculations. The principal new feature in this pair potential 
calculation is the splitting of the lower energy peak into 
bonding and auxtl-bondlng structures. There is no such splitting 
in the higher energy peak, probably because the spread in this 
energy regime caused by the large imaginary part in the surface 

potential related to the inverse of the surface Green function 

' K' makes the two split peaks not resolvable. 

We may embed a cluster of muffin-tin potentials within the 
spherical region I C as shown in figure 6.1 >. The matrix 
representation of the Hamiltonian will then be in terms of a 
basis centred at each of the wells. This is reminiscent of the 
type of approach used in molecular electronic strtjcture 
calculations with an atomic-Uke basis. The calculations will 
then Involve multi-centre integrals. Apart from this 
complication, which has already been dealt with great detail in 
molecular calculations, the rest of the procedure is similar to 


what we have done. 


We have also developed one dimensional analogs of the 
embedding formalism with the lattice structure left Intact. Ve 
have calculated the required Green functions of the effective 
lattice usinc the KKR band theory in one dimension, developed by 
Butler ClS>7d>. Usin^ this we have calculated the local electron 
density of states,embeddln6 a single and a pair muffin-tin 
potential wells in an effective CPA medium. Since these are aU 
model cases we have done the calculations for « -0.3 au and 

V * -0.5 au and c ■ 0.5 and 0.1. 

n 

THis f oi^m^lism is an exact# one. The only appnoximat^ion used 
is trhe Gohenent# Potrent^ial Appnoximat^ion t#Q calculat#e t»he enen^y 
dependent# non-local surface potent^ial of the bulk subsfrafe. 

Let# us now give a brief account# of t#he further application 
of the embedding formalism. The one dimensional embeddini: 
formalism can be generalised to a three dimensional lattice where 
the Green function of the medixim II can be obtained fj>om the 
correspondini^ KKR theory. The difference here is that the three 
dimensional path operator can not be calculated analytically as 
In one dimension. It has to be calculated numerically usin^ 
Ewald^s or other related methods. Computer packages for the f cc 
and bcc lattices are available with us. Once this three 
dimensional formalism is obtained^ calculation of the electronic 
properties of the real materials can be attempted. 


This I'oi'tnallsm can also be applied -Lo study the effects of 
extended defects such as impurity clustering in an alloy and in 
liquid alloys, perhaps short-ranged order. 


Figure 6.1 :A schematic diagram showing i 
of spherically symmetric wells with a 
spherical boundary. 


cluster 
modi f i ed 
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